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We propose a novel method for the effective action of spinor and scalar QED at finite temperature 
in time-dependent electric fields, where charged pairs evolve in a nonadiabatic way. The imaginary 
part of the effective action consists of thermal loops of the Fermi-Dirac or Bose-Einstein distribution 
for the initial thermal ensemble, weighted with factors of the Bogoliubov coefficients for quantum 
effects. And the real part of the effective action is determined by the mean number of produced 
pairs and vacuum polarization at zero temperature. In the weak-field limit, the mean number of 
produced pairs is shown twice the imaginary part. We explicitly find the finite temperature effective 
, action in a constant electric field. 
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I. INTRODUCTION 

In a strong electromagnetic field the vacuum becomes polarized due to the interaction of the electromagnetic field 
^ ' with virtual charged pairs from the Dirac sea. The effective actions in electromagnetic fields have been continuously 
^ investigated since the early work by Sauter, Heisenberg and Euler, and Weisskopf [1] and later on the proper-time inte- 
■ gral for the effective action by Schwinger [2| . The Heisenberg-Euler effective action exhibits both vacuum polarization 
I and pair production, and has many physical applications (for a review and references, see Refs. 0^)- 
^ I The effective actions at zero temperature have been a nontrivial task for general profiles of electromagnetic fields. 
■r^Hj- As a strong electric field always creates pairs from the vacuum, the corresponding effective action contains not only 
, the real part responsible for vacuum polarization but also the imaginary part for the decay of vacuum. Thus, the 
C ■ ' quantum field theory for strong electric fields should properly handle pair creation from the vacuum. The effective 
C5 ' actions have been found for a pulsed-electric field of Sauter type in the resolvent method @ and in the evolution 
\^ • operator method [7,], and the effective action could be found for a spatially localized electric field @. 

However, the QED effective action at finite temperature in electric field backgrounds has been an issue of constant 
: interest and controversy, partly because different formalisms give conflicting results [ol-fl^ and partly because the 
' thermal effects may be important to astrophysical objects involving strong electromagnetic fields. In fact, most 
^ methods for finite temperature field theory may not be directly applied to electric fields due to pair creation from the 
vacuum. The one-loop en ergy momentum tensor of fermions in an initial thermal ensemble was found in a constant 
. , electric field background [l^- Recently the closed-time formalism has been employed to find the QED effective 
' action at finite temperature in 0-f 1 dimension [l6| . The enhancement of pair production by the electric field at finite 
" " ' temperature is also found [l7| . 

The purpose of this paper is two-fold: we first propose a novel method for the effective action at finite temperature 
for time-dependent quantum fields and then find the QED effective action in strong electric field backgrounds. At 
zero temperature the effective action is the scattering amplitude between the out-vacuum and the in-vacuum, which 
is the expectation value of the evolution operator with respect to the in-vacuum Q ■ To extend the in- and out-state 
formalism to finite temperature, we first express the evolution operator in terms of the Bogoliubov coefficients and 
then find the effective action as the expectation value of the evolution operator with respect to the 'thermal vacuum'. 
It turns out that the finite temperature effective action is the trace of the evolution operator weighted with the initial 
thermal ensemble of fermions or bosons, which is equivalent to the 'thermal vacuum' expectation value of the evolution 
operator in thermofield dynamics. The formalism may be applicable to other time-dependent quantum field, whose 
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Hamiltonian H(t) evolves in a nonadiabatic way, that is, out of equilibrium, so that e"'^^^*^ is not the density operator 
satisfying the Liouville-von Neumann equation [ISJ. 

We apply the new method to the QED effective action at finite temperature in time-dependent electric fields. The 
QED effective action consists of the zero-temperature part, the part for thermal and vacuum fluctuations, and the 
finite temperature part without the electric field. The logarithm of the Bogoliubov coefficient plays a role of complex 
chemical potential in the complex thermal distribution for the effective action. The real and the imaginary parts of 
the effective action have an expansion in terms of the Fermi-Dirac or Bose-Einstein distribution and the chemical 
potential. Finally we find the effective action in a constant electric field and discuss it for the Sauter-type electric 
field. 

The organization of this paper is as follows. In Sec. II, we propose a new method to find the finite temperature 
effective action for time-dependent systems. The effective action is given by the trace of the evolution operator and 
the initial density operator, which is equivalent to the expectation value of the evolution operator with respect to the 
thermal vacuum of thermofield dynamics. In Sec. Ill, we find the effective action in spinor and scalar QED in electric 
fields and then elaborate an expansion scheme in terms of the Fermi-Dirac and Bose-Einstein distributions. In Sec. 
IV, we apply the formalism to find the effective action in a constant electric field. Finally, we discuss the controversial 
issue of thermal effects on pair-production rate in Sec. V. 



II. FINITE TEMPERATURE EFFECTIVE ACTION 



We consider both spinor and scalar QED with the time-dependent gauge field A^\{t), which generates a constant 
or time-dependent electric field. For a pulse-like electric field acting for a finite period of time, the ingoing and the 
outgoing vacua are well-defined at Ua = — oo and tout = oo, for which we may choose a gauge A||(iin) = such that 
the ingoing vacuum |0, tin) is nothing but the Minkowski vacuum |0)m- In the case of a constant electric field, we may 
use the asymptotic state as in Ref. The particle and antiparticle have the momentum k and the spin state a, 
whose annihilation operators are denoted by ako-.in and 6kcr,in at t-m = —oo and Oko-.out and 6k<T,out at tout = oo, where 
a = ±1/2 for spinor QED and cr = for scalar QED. Then, the in- and out-vacua are related through the Bogoliubov 
transformations [l^ 

QkCT.out — /^ko-lkcr.in + '^kcr^ko-.in ~ Uka-aka-An^ka-' 

&ko-,out = A'kcr&ko-.in + ^'ko-^^ko-.in — t^kcr^^kcranC^kcr' (1) 

where U^a is the evolution operator, whose form in terms of /iko- and I'ko- is explicitly given in Ref. ^], and the 
coefficients satisfy the relation 

lMk.P + (-l)'+'l'^l|^^k.P = l. (2) 

In the in- and out-state formalism elaborated in Ref. Q, the in- and out-vacua are annihilated by flktr. in/out and 
6kCT, in/out- In fact, the in- and out-vacua are the tensor product of the zero-number states for all k and a. The 
evolution operator transforms the in- vacuum to the out- vacuum as |0,out) = C/|0, in), where U is also the tensor 
product of each C/ko-, that is, U = riko- ^ko-- The zero-temperature effective action per unit volume and per unit time 
is obtained from the scattering amplitude 

e^'d'^dtc^H ^ (o,out|0,in) = (0, in|[/t|0, in). (3) 

Now we extend the zero-temperature effective action to the finite-temperature one for the system with the initial 
density operator. 

Pin = n[^"""^''"''"^l'^k.,in)(nk.,in|], (4) 

ko- 

where /3 = l /fen T, fee being the Boltzmann constant, and i?(k, cr) is the energy for massive charged particles and 
antiparticles. [2^ 

In finite temperature field theories for static systems, one employs either the partition function Z{l3) ~ Tr(/9) or 
the thermal expectation value (0)fl = Tr(Op)/Tr(/9), which is equivalent to the 'thermal vacuum' expectation value, 
(0)/3 = (0, /3, in|O|0, /3, in) |20l - [2a |. However, finite temperature field theories should be modified for time-dependent 
systems since they nonadiabatically evolve the initial states. For such time-dependent quantum fields we propose the 
finite-temperature effective action 
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In fact, the effective action ([5]) is equivalent to (0, /3, in|f/^|0, /3, in) for the 'thermal vacuum' [20l. [2]| 



|0, /3, in) EE -1^ Yl e-"-^^(''--)/ Vka, in) $5 1%., in) 



(6) 



ko- nk<T 



where jfikcr, in) denotes the state for a noninteracting fictitious system of the extended Hilbert space. In fact, Eq. ([5]) 
has the correct zero-temperature limit ([3]) . The effective action ([5|) may be applied to other time-dependent quantum 
fields as well as QED, which evolve in a nonadiabatic way. 



III. QED EFFECTIVE ACTION AT T 

We now advance a method to compute the QED effective action ([5]) in electric fields. Evaluating Eq. ([5]), we obtain 
the effective action at finite temperature per unit volume and per unit time, 

£eff(T, E) ^ (-1)21^1^ J2 [-P^ka + ln(l + (_i)i+2kle-^(-'-^>'-)) _ ln(I + (_i)i+2kle-/3-k)j ^ (7) 



where Wk = ym? -|- + (fc|| + qA\\)^ and 

^=6^"''% {Zka = zr{k,a)+iz,{k,a)). (8) 

The summation is over all possible states such as momenta and spin states. Each term in Eq. ([T]) has the following 
interpretation: the first term is the effective action Ccs{T — 0, E) at zero temperature, the second term is the combined 
effect of thermal and quantum fluctuations, while the last term is the subtraction of the effective action (potential 
energy) Ccs{T^ E — Q) at finite temperature without the electric field. From now on we subtract the zero-temperature 
part from the effective action and let 

^C,s{T,E)^C,s{T.E)~C,fi{Q,E). (9) 

Note that Zkcr(i?), which depends on the electric field E and Zkcr(O) = 0, plays a role of complex chemical potential, 
as will be explained below. 

Further, we elaborate an expansion scheme for the effective action in terms of the Fermi-Dirac or Bose-Einstein 
distribution and Zko-. First, the imaginary part of the effective action ([9]) can be expanded as 

Im(A/:eff) = (-I)i+2|'^li5]f;Ml^^^^Z^[(e^- - I)^- + (e^^- - 1)^], (10) 

kcr J = l ^ 

where rtF/B(k) denotes cither the Fermi-Dirac distribution ni?(k) = I/(e'^"'' -I- 1) for spinor QED or nB(k) — I/(e^"'' — 
1) for scalar QED. Second, the real part of the effective action © is given by 

Re(A/:eff) = (-l)'l'^l E E ^ ■■ sinO-/3^.(k, a)). (11) 

kcr J = l ^ 

The thermal factors in Eqs. (jlOp and (jlip correspond to thermal loops in the diagrammatic representation, which are 
weighted with factors from quantum fluctuations 

Now, we give physical interpretations for the effective action. In the weak-field limit {qE <^ m?) where /3zi(k) <C 1, 
the real part (fTTj) approximately is 



/3zi(k,a-) 

ko- 

while the imaginary part (|10p approximately leads to 

2Im(A/:eff) « (-l)'l"l W<r?np/B{\i). (13) 

kcr 
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Thus, the imaginary part may be regarded as the pair-production rate due to thermal and quantum effects. The 
thermal effects suppress the fermion pair production due to the Pauli blocking but enhance the boson pair production 
due to the Bose-Einstein condensation, as expected. In Ref. [l^ , the mean number of produced pairs with a given mo- 
mentum k at T is given by A/'"P(r) = Y.\,a l^k^J^ tanh(^Wk/2) for spinor QED and A/'"'=(T) = Y.k^ l^^icrP coth(/3wk/2) 
for scalar QED. So the mean number, AM — {Af{T) —Af{0))/2, of one species of particle or antiparticle due to thermal 
effects approximately satisfies the relation between the mean number and the imaginary part: 

AA^ = ^ Wwafnp/BO^) « 2Im(A£eff). (14) 

The relation between the mean number of produced pairs and twice of the imaginary part also holds at T = in the 
weak-field limit 

A few comments are in order. The series of the real part (fTTj) may be summed as [l^ 

r sin(/3zi(k)) ] 

Re(A/:eff) = garctan[ ^^^^^_^^^,^^ ^ ^^^^^^^^^^^^ j ■ (15) 

Using £cff(0, E) — (— 1)^^^''^'* X^ka /'■^ko- from Eq. ([7]) and the Bogoliubov relation ([2]), we have the real and imaginary 
parts 

/3z.(k,a) = (-l)i+2klini(£^^(0,£;)) = -i±|Min(l + (_l)2kl|^^^|2)^ 

/3z,(k,a) = (-l)2|-lRe(/:eff(0,i;)). (16) 

Then the effective action (jlSp at finite temperature can be written in terms of the mean number, the vacuum 
polarization, and the thermal distribution as 

Re(A£.«.) = (-1)21-1 Varctanf sin(Re(£.40, i^))) 1^ 

Lg/3.,(i + (_^)2k||^^^|2)^ + (_l)i+2k| cos(Re(/:cff(0,£;)))J 

Another interesting observation is that the mid-term in Eq. 

WasiT, E) = (-l)i+2|'^lz ^ ln(l + (_i)i+2kle-/3("k-3k.))^ (18) 

is reminiscent of the potential energy and carries both thermal and quantum effects. Equation suggests Zko- 
as the chemical potential and the variation with respect to z-^„ yields the Fermi-Dirac or Bose-Einstein distribution. 



IV. APPLICATIONS 



In this section we find the QED effective action in a constant electric field and discuss a Sauter-type electric field, 
E{t) — Easech^ {t / t) in Ref. [7J. In the time-dependent gauge, A\\^(t) = —Et for the constant electric field and 
A\\(t) — -~Eot{1 + tanh(i/T)) for the Sauter-type electric field, the energy of charged particles explicitly depend on 
time. We will take the weak-field limit {qE <C m?), where the real part and the imaginary part of the 
approximate effective action can be worked out for the constant electric field and in principle for the Sauter-type 
electric field. 

In the constant electric field, the state along the direction of the electric field is asymptotically determined, whose 
momentum integral gives a factor qE/{2TT) 0. Using the mean number of produced pairs. 



which is independent of the spin states, the imaginary part (|13p is given by 



1 



3m{qE) 



n=0 
2 



(19) 



J7j2g-/3m(n-|-l) 

27rm2 -|- l3mqE{n + 1) 
3l3m{qEf 

(27rm2 -I- /3mqE{n + 1))2 {2Trm? + f3mqE{n + 1))3 



(20) 
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The factor in front of the summation is the leading term of the imaginary part at zero temperature. Further, in the 
low-temperature limit (/3m 3> 1), the leading term of Eq. (I^Hl) is 



\2ttJ 



m2 e"^™ 



qE /3m 



qE 



(21) 



In the special case of thermal effect dominance, neglecting all terms of m/PqE, the first series in Eq. ((20|) approximately 
leads to 



Similarly, using Re(£off (0, i?)) in Ref. j?:], the real part p^ . for instance, of spinor QED is given by 



(22) 



n=0 



(=2 



(20! V27r/ m4'-2 



3 - 21, a) + ^"^(" + ^^^(3, 5 - 21, a) - ^-^^^{4, 6 - 21, a) 
V 4 8 



(*(1, 3 ~ 21, 7) + ^^(3, 5 ~ 2/, 7) - ^^(4, 6 - 2/, 7) + •••)+ ' 



where B21 is the Bernoulli number, ^1/ denotes the second confiuent hypergeometric function [25j . and 

/3m{n + 1) f3m{n + 1) 7rm^ 

a = z , 7 = 



2 ' ' 2 g£: 

In the low-temperature limit (/3m ^ 1), the series of / = 2 in Eq. (1^^ leads to 



ReiACZ{T,E)) 



(2^/g^y 
' 45m4 V 27r / 



(-1)— ln(l + e-^") 
/3m 



(23) 



(24) 



(25) 



Here the factor in front of the square bracket is the real part at zero temperature. The real part of effective action in 
scalar QED may be found in a similar way. 

Finally, we discuss the Sauter-type electric field. The charged particle has the free energy Wk.in — V rn^ + before 
the onset of the electric field while it has Wk.out — \/rn^ + -I- [k^ — 2qEoT)'^ after the completion of the interaction. 
At zero temperature, the mean number of produced pairs, Eqs. (68) and (83), and the vacuum polarization, Eqs. 
(66) and (80) of Ref. [7|, which depend on Wk.inj Wk,out, and A = yj {qE^r'^Y ~ {2\cj\ — 1)^/4, lead to the effective 
action, Eqs. ((T^ and (IT5|) . To find analytical expressions for the effective action would be more complicated than the 
constant electric field, which will be addressed elsewhere. 



V. CONCLUSION 



In this paper we have advanced a new method for the finite-temperature effective action for time-dependent quantum 
fields and have studied the one-loop effective action of spinor and scalar QED at finite temperature in a constant or 
time-dependent electric fields. External electric fields make the vacuum unstable against pair production, which is 
a consequence of the out- vacuum differing from the in- vacuum. The instability enforces a careful application of the 
finite-temperature field theory to time-dependent quantum fields. The finite-temperature effective action ([5]) is given 
by the trace of the initial thermal ensemble evolved by the time-evolution operator, which is equivalent to the thermal 
vacuum expectation value of the evolution operator in thermofield dynamics. 

The imaginary part (jlOp of the effective action exhibits factorization into thermal factors and quantum factors, 
which correspond to thermal loops in the diagrammatic representation with vertices of the external electric field. In 
the weak-field limit [qE <^ m?), twice of the imaginary part is the mean number of produced pairs, as shown in Eq. 
([T3l) . However, the thermal and quantum effects are intertwined in the real part of the effective action, Eqs. (|11|) . 
([T5|). and P7|) . In fact, the finite-temperature effective action fTT]) is determined by the vacuum polarization at zero 
temperature, the mean number of produced pairs, and thermal distribution. In the weak-field and lower-temperature 
limits, the leading factors of the real and imaginary parts, (12^ and (I^Sj) . of the effective action in a constant electric 
field are proportional to those at zero temperature and the potential energy for the rest mass in spinor and scalar 
QED. 
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Our results show many interesting aspects. First, the imaginary part (jlOp of the thermal contribution does not 
vanish for any non-zero electric field, which implies thermal effects on pair production and thus may resolve the 
controversial issue of thermal effects on pair production: thermal effects are shown to exist in Refs. [13, while 
in Refs. jTlJ no thermal effects are found. Further, in the weak-field limit for small pair production, twice of the 
imaginary parts (|14p are the pair-production rate at T, which was shown in Rcf. (loj . Though our in- and out-state 
formalism differs from the imaginary-time formalism, the imaginary part ()13p in the weak-field limit is the pair- 
production rate times the Fermi-Dirac or Bose-Einstein distribution, which may correspond to two-loop dominance 
in Ref. Second, the Bogoliubov coefficient ([5]), which is responsible for vacuum polarization at T = 0, plays a 

role of chemical potential in the effective action ([7]) and in the potential energy (ITSl) at T. In fact, the variation of 
the effective action with respect to the chemical potential yields the Fermi-Dirac or the Bose-Einstein distribution. 
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